Naval  Research  Laboratory 

Stennis  Space  Center,  MS  39529-5004 


NRL/FR/7322-96-9656 


Tidal  Elevations  in  the  Yellow  and  East 
China  Seas  From  a  Data-Assimilative, 
Barotropic  Tidal  Model 


Cheryl  Ann  Blain 

Ocean  Dynamics  and  Prediction  Branch 
Oceanography  Division 


July  3,  1997 


19970128  179 


Approved  for  public  release;  distribution  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OBMNo.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  or  any  other  aspect  of  this  collection 
of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Sendees,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. _ 


1 .  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

July  3,  1997  Final 


2.  REPORT  DATE 

July  3,  1997 


4.  TITLE  AND  SUBTITLE 


Tidal  Elevations  in  the  Yellow  and  East  China  Seas  From  a  Data-Assimilative, 
Barotropic  Tidal  Model 


6.  AUTHOR(S) 


5.  FUNDING  NUMBERS 


Job  Order  No. 


573509706 


Program  Element  No.  0603207N 
Project  No.  X2008 


Cheryl  Ann  Blain 


Accession  No. 


DN252022 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

Oceanography  Division 

Stennis  Space  Center,  MS  39529-5004 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NRL/FR/7322-96-9656 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Space  and  Naval  Warfare  Systems  Center 

PMW-165 

Washington,  DC 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

The  Yellow  Sea  is  a  very  shallow-water  body  with  a  strong  tidal  signal.  Its  strategic  importance  had  led  the  U.S.  Navy  to 
study  tidal  response  in  the  Yellow  and  East  China  Seas  using  the  Colorado  University  (CU)  data-assimilative,  barotropic  tidal 
model.  The  CU  model  has  been  adapted  to  accept  open-ocean  boundary  forcing  from  the  Grenoble  global  tidal  data  base,  which 
is  considered  more  reliable  than  tides  currently  obtained  from  the  CU  global  tidal  model.  Furthermore,  the  number  and  quality 
of  the  data  assimilation  stations  utilized  by  the  CU  Yellow  Sea  model  is  improved.  A  detailed,  semi-automated  procedure  for 
assigning  data  stations  for  assimilation  into  the  CU  model  domain  has  been  developed. 

Characterization  of  the  tidal  elevations  in  the  Yellow  and  East  China  Seas  from  CU  model  simulations  demonstrate  that  the 
Mj  constituent  is  the  dominant  tidal  frequency  in  the  Yellow  Sea  with  amplitudes  of  nearly  2  m  along  the  coast  of  west  Korea. 
Ki  is  the  most  significant  diurnal  component  with  amplitudes  1/4  of  the  Mj  tide. 

A  series  of  experiments  is  conducted  to  test  the  influence  of  the  location  of  assimilated  data  on  predictions  of  the  dominant 
tidal  constituents  Mj  and  K,.  From  an  analysis  of  four  data  station  groupings,  coastal,  shallow,  deep,  and  open  boundary  located 
stations,  observations  located  on  the  continental  shelf  away  from  the  coastline  are  most  influential  in  the  prediction  of  primary 
tides,  Mj  and  Kj  in  the  Yellow  Sea.  Coastal  station  data  may  be  more  influential  in  the  prediction  of  nonlinear  tidal  interactions, 
since  these  gauges  are  often  located  in  harbors  or  inlets  where  nonlinear  tidal  interactions  are  prominent  and  wave  action  may 
be  present.  Furthermore,  coastal  data  is  assimilated  into  the  model  at  the  nearest  nodal  locations  introducing  an  additional  source 
of  error.  This  examination  of  the  data  assimilation  component  of  the  CU  model  should  enhtmee  development  of  data  sampling 
strategies  to  improve  future  tidal  model  predictions. 


14.  SUBJECT  TERMS 


Yellow  Sea,  shallow-water  modeling,  tidal  modeling 


15.  NUMBER  OF  PAGES 

24 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

Unclassified 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 


20.  LIMITATION  OF  ABSTRACT 


Same  as  report 


NSN  7540-01-280-5500 


1 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std.  Z39-18 
298*102 
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but  phases  are  poorly  represented.  RMS  errors  for  the  data-assimilative  simulation  are  lower  than  those  for  the  nonassimilative 
run,  but  comparison  of  absolute  tidal  amplitude  errors  indicate  overdamping  by  the  CU  model. 
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EXECUTIVE  SUMMARY 


The  Yellow  Sea  is  a  very  shallow-water  body  with  a  strong  tidal  signal.  Its  strategic  importance 
had  led  the  U.S.  Navy  to  study  tidal  response  in  the  Yellow  and  East  China  Seas  using  the  Colorado 
University  (CU)  data-assimilative,  barotropic  tidal  model.  The  CU  model  has  been  adapted  to 
accept  open-ocean  boundary  forcing  from  the  Grenoble  global  tidal  data  base,  which  is  considered 
more  reliable  than  tides  currently  obtained  from  the  CU  global  tidal  model.  Furthermore,  the  num¬ 
ber  and  quality  of  the  data  assimilation  stations  utilized  by  the  CU  Yellow  Sea  model  is  improved. 
A  detailed,  semi-automated  procedure  for  assigning  data  stations  for  assimilation  into  the  CU 
model  domain  has  been  developed. 

Characterization  of  the  tidal  elevations  in  the  Yellow  and  East  China  Seas  from  CU  model 
simulations  demonstrate  that  the  M2  constituent  is  the  dominant  tidal  frequency  in  the  Yellow  Sea 
with  amplitudes  of  nearly  2  m  along  the  coast  of  west  Korea.  Ki  is  the  most  significant  diurnal 
component  with  amplitudes  1/4  of  the  M2  tide. 

A  series  of  experiments  is  conducted  to  test  the  influence  of  the  location  of  assimilated  data 
on  predictions  of  the  dominant  tidal  constituents  M2  and  Ki.  From  an  analysis  of  four  data  station 
groupings,  coastal,  shallow,  deep,  and  open  boundary  located  stations,  observations  located  on  the 
continental  shelf  away  from  the  coastline  are  most  influential  in  the  prediction  of  primary  tides,  M2 
and  Ki  in  the  Yellow  Sea.  Coastal  station  data  may  be  more  influential  in  the  prediction  of  non¬ 
linear  tidal  interactions,  since  these  gauges  are  often  located  in  harbors  or  inlets  where  nonlinear 
tidal  interactions  are  prominent  and  wave  action  may  be  present.  Furthermore,  coastal  data  is 
assimilated  into  the  model  at  the  nearest  nodal  locations  introducing  an  additional  source  of  error. 
This  examination  of  the  data  assimilation  component  of  the  CU  model  should  enhance  development 
of  data  sampling  strategies  to  improve  future  tidal  model  predictions. 

Validation  of  the  CU  tidal  model  is  accomplished  by  comparing  predicted  and  observed  tidal 
elevations  at  49  locations.  For  eight  primary  tidal  constituents,  absolute  and  root-mean- square 
(RMS)  amplitude  and  phase  errors  are  examined  for  data-assimilative  and  nonassimilative 
computations.  CU  model  results  for  the  M2  and  Ki  tides  are  within  13%  of  measured  values,  but 
phases  are  poorly  represented.  RMS  errors  for  the  data-assimilative  simulation  are  lower  than  those 
for  the  nonassimilative  run,  but  comparison  of  absolute  tidal  amplitude  errors  indicate  overdamping 
by  the  CU  model. 

A  look  at  the  sensitivity  of  the  bottom  friction  coefficient  specification  on  model  computations 
shows  that  the  M2  tide  is  most  affected.  An  increasing  friction  coefficient  decreases  tidal  amplitudes, 
although  absolute  errors  indicate  an  optimal  friction  factor  between  0.002  and  0.003,  which  is  in 
good  agreement  with  published  values  for  the  Yellow  Sea. 
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TIDAL  ELEVATIONS  IN  THE  YELLOW  AND  EAST  CHINA  SEAS  FROM  A 
DATA-ASSIMILATIVE,  BAROTROPIC  TIDAL  MODEL 


1.0  INTRODUCTION 

Increased  interest  by  the  U.S.  Navy  in  the  prediction  of  tides  in  coastal  areas  and  semi-enclosed 
seas  is  a  direct  result  of  a  shift  in  emphasis  from  deep-water,  antisubmarine  warfare  to  littoral 
warfare  scenarios.  In  recent  times,  the  Yellow  Sea  has  grown  in  strategic  importance  due  to  the 
dynamic  political  situation  between  North  and  South  Korea  and  China  and  Taiwan,  countries  whose 
borders  all  ring  the  Yellow  and  East  China  Seas.  As  a  consequence,  the  tidal  response  in  the  Yellow 
and  East  China  Seas  is  of  interest  to  military  planners  in  the  U.  S.  Navy.  The  Yellow  Sea  itself  is 
quite  shallow  and  its  circulation  is  largely  dominated  by  a  strong  tidal  signal.  Accurate  prediction 
of  the  tidal  elevation  is  essential  to  the  success  of  any  military  operation  in  this  region. 

Implementation  of  a  numerical,  barotropic  tidal  model  for  any  water  body  involves  specification 
of  bathymetry,  boundary  forcing,  and  model  parameters  such  as  the  bottom  friction  coefficient. 
Each  of  the  data  sources  for  these  model  components  has  associated  uncertainties,  particularly  in 
coastal  areas  and  marginal  seas.  Bathymetry  is  commonly  known  only  to  1/12°  and  can  be  inaccu¬ 
rate  in  near-shore  regions.  Boundary  forcings  are  taken  from  sparse  data  or  sea-surface  heights, 
which  themselves  have  been  computed  from  larger  scale  ocean  models.  Model  parameter  selection 
is  often  based  on  experience  or  values  are  “tuned”  during  model  calibration. 

To  partially  compensate  for  uncertainties  in  model  data  and  overcome  errors  due  to  inadequate 
grid  resolution,  observed  tidal  heights  are  merged  with  numerical  tidal  models  through  data  assimilation. 
In  coastal  zones  and  marginal  seas  where  complex  bathymetry  and  irregular  shorelines  produce 
highly  nonlinear  tidal  interactions,  the  assimilation  of  tidal  heights  into  model  computations  may 
significantly  improve  model  performance.  Furthermore,  for  short-term  forecasts,  the  assimilation  of 
observed  sea-surface  heights  on  real-time  scales  may  be  required  to  achieve  a  desired  accuracy  in 
computed  tidal  elevations. 

The  U.  S.  Navy  currently  exercises  a  relocatable,  data-assimilative,  barotropic  tidal  model  in 
the  Yellow  Sea.  This  model  is  the  Colorado  University  (CU)  tidal  model  and  is  a  variant  of  the 
Princeton  finite  difference,  hydrodynamic  model.  The  CU  model  utilizes  a  coarse  grid  and  relies 
heavily  on  data  assimilation  to  achieve  accurate  tidal  predictions  (Kantha  et  al.  1994).  In  the 
Yellow  and  East  China  Seas,  ample  observation  points  are  available  for  assimilation  into  the  tidal 
model  or  for  comparison  during  model  evaluation.  Data  sources  include  tide  gauges  as  well  as 
TOPEX/Poseidon  altimeter  measurements  (Lee-Lueng  et  al.  1994). 

From  CU  model  predictions  of  tidal  elevation,  a  spatial  description  of  tidal  forcing  in  the 
Yellow  and  East  China  Seas  is  obtained.  Cotidal  charts  for  eight  primary  tidal  constituents  are 
derived  and  their  relative  contributions  to  the  overall  tidal  signal  are  assessed.  The  influence  of  the 
location  of  assimilated  data  on  computed  sea-surface  height  predictions  for  the  major  semidiurnal, 
M2,  and  diurnal,  Ki,  tidal  constituents  is  demonstrated.  Tidal  data  located  at  114  stations  are 
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divided  into  four  groups  corresponding  to  the  bathymetry  at  their  location  or  their  proximity  to 
open-ocean  boundaries.  Tidal  heights  computed  by  assimilating  each  group  of  data  stations  separately 
are  qualitatively  compared  to  assess  how  differences  in  the  tidal  solutions  can  be  attributed  to 
station  location.  Finally,  the  CU  model  capability,  with  respect  to  prediction  of  tidally  induced  sea- 
surface  height,  is  quantitatively  evaluated.  The  1 14  locations  of  measured  tides  are  separated  into 
assimilation  stations  and  observational  data  locations.  Absolute  and  root-mean-square  (RMS)  errors 
of  model  computations  derived  from  comparisons  to  the  observed  tidal  elevations  are  presented.  As 
part  of  this  model  evaluation,  the  sensitivity  of  bottom  friction  coefficient  on  model  results  is 
briefly  considered. 


2.0  NUMERICAL  MODELING  STRATEGY 

The  hydrodynamic  computations  are  made  using  a  fully  nonlinear,  two-dimensional,  barotropic 
tidal  model  developed  at  CU  (Kantha  et  al.  1994;  Kantha  1995).  The  governing  equations  are 
derived  from  the  three-dimensional  equations  for  mass  and  momentum  balance  subject  to  the  hydrostatic 
assumption,  the  Boussinesq  approximation,  and  vertical  integration  over  the  water  column.  Written 
in  spherical  coordinates  the  equations  take  the  following  form: 
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R  cos(l) 


a  _  a  _ 

^RDux  +^Rcos<?Dui^ 


a  ^  tan  ()) 

■^^Dux+  u-Vux -fDu^  + - - - : 


(2) 


+  u-Vl^ -fDui - - - =  (3) 

where  t  represents  time.  A,,  ({)  are  degrees  longitude  (east  of  Greenwich  positive)  and  degrees  latitude 
(north  of  the  equator  positive),  T|  is  the  free  surface  elevation  relative  to  the  geoid,  ux  >  ^  the 
depth-averaged  horizontal  velocities,  R  is  the  radius  of  the  Earth,  D  is  the  bathymetric  depth 
relative  to  the  geoid,  /  =  2Qsin  (|)  is  the  Coriolis  parameter,  £2  is  the  angular  speed  of  the  Earth, 
p  is  the  atmospheric  pressure  at  the  free  surface,  g  is  the  acceleration  due  to  gravity,  P  is  the 
reference  density  of  water,  is  the  effective  Newtonian  equilibrium  tidal  potential,  x^x,  Xbt^  are 
the  bottom  stresses,  and  Fx,  are  the  horizontal  mixing  terms. 

In  the  current  implementation  of  the  CU  model,  both  baroclinic  forcings  and  wind  stress  are 
neglected.  Bottom  friction  terms  are  parameterized  using  a  quadratic  drag  law  and  a  Laplacian 
formulation  characterizes  the  horizontal  mixing.  The  viscosity  coefficient  in  the  horizontal  mixing 
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terms  is  composed  of  two  parts:  a  constant  and  a  Smagorinsky  nonlinear  portion,  which  is  dependent 
on  mean  strain  rate  (Smagorinsky  1963). 

The  equilibrium  tidal  potential,  which  represents  astronomical  forcing  by  celestial  bodies, 
the  sun  and  the  moon,  is  taken  from  Schwiderski  (1980): 

^  (?L,(p.^)  =  ^  ^jn h  (9)  COS  {IsttlTjn  +jX  +  ,  (4) 

nij 

where  Bjn  is  the  amplitude  of  constituent  n  for  species  j  =  0,  1,2,  the  corresponding  period  is  Tjn, 
Xjn  is  the  astronomical  argument,  and  Lj  (tp)  is  defined  by: 


Bq  (9)  =  3  sin^tp  -  1, 

(5) 

Li  (9)  =  sin29. 

(6) 

L2  (9)  =  cos,\ 

(7) 

The  use  of  frequency  and  time-dependent  nodal  factors  is  not  required  by  the  CU  model  when 
using  the  following  methodology.  Primary  constituents,  four  semidiurnal  (M2,  S2,  N2,  and  K2),  four 
diurnal  (Ki,  Oi,  Pi,  and  Qi),  and  three  long-term  (Mf,  Mm,  Ssa),  are  simulated  individually  using 
the  CU  model.  The  orthotide  approach  of  Groves  and  Reynolds  (1975)  can  then  be  applied  to 
compute  the  total  tidal  response,  which  includes  60  high-frequency  components  and  15  of  the  most 
significant  long-term  equilibrium  tides.  The  only  time-dependent  information  needed  for  a  model 
simulation  is  the  date  and  year.  Basis  for  the  orthotide  approach  lies  in  the  representation  of  tides 
as  a  linear  combination  of  a  special  set  of  orthogonal  tidal  functions.  The  reader  is  referred  to 
Groves  and  Reynolds  (1975)  or  Kantha  (1995)  for  a  more  complete  discussion  of  the  implementation 
and  facility  of  the  orthotide  method. 

The  governing  equations  are  discretized  using  an  explicit,  finite  difference  scheme  over  an 
orthogonal,  curvilinear  grid.  The  grid  is  known  as  the  staggered  Arakawa-C  grid  where  sea-surface 
height  and  velocity  are  computed  at  locations  within  the  grid  cell  that  differ  by  one  half  the  grid 
spacing.  In  time,  a  forward  differencing  scheme  is  implemented  to  reduce  stability  restrictions  on 
the  time  step.  Numerical  formulations  in  the  CU  tidal  model  are  nearly  identical  to  that  of  the 
well-known  Princeton  hydrodynamic  model  (e.g.,  Blumberg  and  Mellor  1987;  O’Connor  1991). 

A  distinctive  feature  of  the  CU  model  over  other  shallow-water  models  (e.g.,  Blumberg  and 
Mellor  1987;  Lynch  and  Werner  1991;  Luettich  et  al.  1992)  is  the  assimilation  of  tidal  component 
data  taken  from  coastal  tide  gauges  and/or  altimetry  sources.  A  rather  simplistic  data-assimilation 
procedure,  nudging,  is  implemented.  Nudging  is  a  localized  version  of  optimal  interpolation  tech¬ 
niques  that  are  widely  used  in  conjunction  with  dynamical  models  having  heterogeneous  observations 
(Ghil  and  Malanotte-Rizzoli  1991).  Model  predicted  sea-surface  heights  at  assimilation  locations 
are  replaced  at  predetermined  intervals  by  a  weighted  sum  of  the  model  prediction  and  the  observed 
sea-surface  height.  Weights  for  the  data  are  determined  a  priori  and  are  selected  to  reflect  a  level 
of  confidence  in  the  accuracy  of  assimilated  observations. 

2.1  Upgrades 

The  CU  model  has  been  adapted  to  accept  open-ocean  boundary  forcing  from  the  Grenoble 
global  tidal  data  base  (LeProvost  et  al.  1994),  which  is  considered  more  reliable  than  tides  currently 
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implemented  and  obtained  from  the  CU  global  tidal  model  (Kantha  1995).  The  Grenoble  tide 
solutions  are  derived  from  a  finite  element  hydrodynamic  model.  For  the  M2  and  S2  tidal  constitu¬ 
ents,  empirically  computed  tidal  heights  derived  from  TOPEX/Poseidon  data  are  assimilated  into 
the  model  computations  for  an  improved  solution.  The  Grenoble  tidal  solutions  have  a  resolution 
down  to  10  km  along  the  coasts  and  the  released  data  base  is  gridded  with  a  0.5-deg  resolution  for 
convenience. 

Using  a  constructed  preprocessor,  any  or  all  of  eight  primary  tidal  constituents  are  directly  accessed 
from  the  Grenoble  tidal  data  base  and  then  interpolated  onto  open-ocean  boundary  points  associated 
with  the  CU  model  grid.  Tidal  constituent  data  along  the  open-ocean  boundaries  takes  the  form  of 
amplitudes  and  phases  and  is  written  in  a  format  that  interfaces  with  the  CU  model  setup. 

Another  improvement  to  the  CU  model  is  the  inclusion  of  more,  high-quality  data-assimilation 
stations  in  and  around  the  Yellow  Sea.  A  total  of  114  locations  throughout  the  Yellow  and  East 
China  Seas  are  identified  as  either  having  tide  gauge  observations  or  TOPEX/Poseidon-derived  sea- 
surface  heights.  For  quality  control  purposes,  the  tide  gauges  included  in  this  group  were  active  for 
at  least  180  days  and  are  extracted  from  the  International  Hydrographic  Office  data  (IHO  1979). 
Altimetric  tides  are  derived  from  TOPEX  data  according  to  Desai  and  Wahr  (1994). 

A  detailed,  semi-automated  procedure  is  created  to  assign  data  station  locations  for  assimilation 
into  the  CU  model  domain.  Given  an  assimilation  station  data  file  in  the  form  specified  by  the 
CU  model,  domain  coordinates  in  latitude  and  longitude,  grid  spacing,  and  bathymetry,  data  stations 
are  located  with  the  domain  or  flagged  if  they  are  not  in  the  domain  or  are  located  over  land. 
Flagged  stations  can  be  manually  included  as  assimilation  stations  by  specifying  a  nearby  nodal 
location.  Other  stations,  located  over  water  and  within  the  domain,  are  automatically  assigned  to 
the  nearest  grid  point  for  assimilation. 


3.0  YELLOW  SEA  TIDAL  MODEL 

The  Yellow  Sea,  depicted  in  Fig.  1,  is  a  very  shallow,  semi-enclosed  sea  bordered  by  China  to 
the  north  and  west  and  Korea  in  the  east.  The  southern  Yellow  Sea  merges  into  what  is  known 
as  the  East  China  Sea.  As  seen  in  Fig.  1,  bathymetry  over  much  of  the  Yellow  Sea  is  less  than  80  m. 
The  200-m  isobath  in  the  East  China  Sea  forms  a  sill  beyond  which  depths  drop  sharply  to 
1000  m.  The  source  of  bathymetry  is  the  digital  bathymetric  data  base  (DBDB-5)  on  which  a 
minimum  depth  of  5  m  is  imposed  to  eliminate  drying  of  computational  points  along  the  shoreline. 
The  domain  selected  for  these  experiments  includes  both  the  Yellow  and  the  East  China  Seas, 
extending  as  far  south  as  Taiwan  and  eastward  to  the  Japanese  island,  Kyushu.  The  finite  difference 
computational  grid  for  this  domain  is  rectangular  with  a  relatively  coarse,  regular  spacing  of  0.2  deg 
(~20  km).  Of  the  6106  computational  points,  3557  are  over  water  and  considered  active. 

Eight  primary  tidal  constituents,  four  semidiurnal  (M2,  S2,  N2,  K2)  and  four  diurnal  (Ki,  Oi, 
Pb  Ql)  tides,  are  the  only  forcing  considered  in  the  Yellow  Sea.  Recall  the  modeling  strategy  is 
such  that  each  constituent  is  simulated  separately  within  the  tide  model.  The  eight  constituents  are 
forced  independently  on  the  interior  of  the  domain  through  the  tidal  potential  term.  At  the  southern 
and  eastern  boundaries,  values  for  the  tidal  constituents  are  prescribed  using  results  from  the 
Grenoble  global  tide  model  (LeProvost  et  al.  1994).  Tidal  heights  along  the  Tsushima  Strait  at 
the  entrance  to  the  Sea  of  Japan  are  taken  from  Kantha  (1994). 


Model  parameters  are  identical  for  all  experiments  unless  specified  otherwise  so  that  comparisons 
between  model  simulations  is  possible.  The  bottom  friction  coefficient  is  uniform  and  equal  to 
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0.002  for  all  model  runs  unless  otherwise  indicated.  The  constant  part  of  the  horizontal  viscosity 
coefficient  has  a  value  of  0.02  and  the  nonlinear  portion  associated  with  the  Smagorinsky  formu¬ 
lation  is  specified  as  100  m^/s.  Simulations  are  spun  up  from  a  homogeneous  initial  condition  using 
a  ramp  in  time  equivalent  to  five  inertial  periods.  Model  computations  proceed  for  a  total  of  10  days 
using  a  time  step  of  24  s.  Long-term  simulations  of  150  days  produce  cotidal  charts  identical  to 
those  computed  using  10  days  of  computation. 

A  total  of  114  locations  throughout  the  Yellow  and  East  China  Seas  are  identified  as  either 
having  tide  gauge  observations  or  TOPEX/Poseidon-derived  sea-surface  heights.  For  quality  control 
purposes,  the  tide  gauges  included  in  this  group  were  active  for  at  least  180  days  and  are  extracted 
from  the  International  Hydrographic  Office  (IHO)  data  (IHO  1979).  Altimetric  tides  are  derived  from 
TOPEX  data  according  to  Desai  and  Wahr  (1994). 


6 


Cheryl  Ann  Blain 


4.0  TIDAL  ELEVATIONS  IN  THE  YELLOW  SEA 

Two  CU  tidal  model  simulations  compute  the  amplitudes  and  phases  of  eight  major  tidal 
constituents,  M2,  S2,  N2,  K2  and  Ki,  Oi,  Pi,  Qi,  in  the  Yellow  and  East  China  Seas.  The  first, 
denoted  Run  A,  includes  the  assimilation  of  observed  tidal  heights  at  65  stations,  whereas  the 
second.  Run  NA,  includes  no  data  assimilation.  In  general,  tides  in  the  Yellow  Sea  are  influenced 
by  lateral  boundaries  and  depths  less  than  200  m.  Tidal  waves  propagate  from  the  deep  sea  onto 
the  continental  shelf  where  they  are  amplified  along  the  west  Korean  coast  and  distorted  upon 
entering  the  Bohai  Sea.  Tidal  amplitudes  in  the  East  China  Sea  tend  to  decrease  in  a  direction 
towards  the  Sea  of  Japan. 

Mean  amplitudes,  calculated  over  all  active  node  points,  of  all  eight  major  tidal  components 
from  model  simulations  A  and  NA,  are  presented  in  Table  1.  The  semidiurnal  M2  tide  dominates 
the  tidal  response  in  the  Yellow  Sea.  The  second-largest  contribution  to  the  tidal  signal  comes 
from  the  S2  tidal  frequency,  which  has  magnitudes  ranging  from  39-46%  of  the  M2  tidal  amplitude. 
Amplitudes  of  the  strongest  diurnal  component,  Ki,  are  slightly  greater  than  1/4  of  those  associated 
with  the  M2  constituent.  The  five  remaining  semidiurnal  and  diurnal  constituents  comprise  a 
comparatively  small  part  of  the  overall  tidal  amplitude  in  the  Yellow  Sea. 

When  no  data  is  assimilated  (Run  NA),  model-predicted  tidal  heights  are  generally  larger  as 
reflected  in  the  greater  mean  amplitude  when  compared  to  that  of  Run  A.  However,  for  both 
simulations  NA  and  A,  the  relative  magnitudes  between  each  of  the  eight  tidal  constituents  remain 
stable  with  only  some  variation  in  the  S2  and  N2  magnitudes. 

Cotidal  charts  for  the  dominant  semidiurnal  and  diurnal  tides,  M2  and  Ki,  are  shown  in  Figs.  2 
and  3.  While  not  shown,  cotidal  charts  for  the  remaining  six  semidiurnal  and  diurnal  tides  exhibit 
similar  patterns  as  those  for  M2  and  Ki,  respectively.  These  charts,  derived  from  Run  A  model 
computations,  include  the  assimilation  of  tidal  data  at  65  locations  (shown  later  in  Fig.  7).  The  M2 
frequency  has  four  amphidrome  locations  or  points  of  zero  amplitude  from  which  tidal  phase  lines 
radiate  as  shown  by  the  dotted  lines  in  Fig.  2.  Extreme  values  of  the  M2  tide  are  nearly  200  cm 
along  the  west  Korean  coast. 

For  the  Ki  tidal  component,  two  amphidromes  are  evident  in  dotted  phase  lines  plotted  in 
Fig.  3.  Again,  tidal  heights  for  Ki  are  slightly  greater  than  1/4  of  those  associated  with  M2  and 
maximum  amplitudes  of  nearly  40  cm  are  computed  on  the  west  coast  of  Korea  and  in  the  Bohai 
Sea.  These  regions  of  tidal  amplification  are  similar  to  those  seen  in  the  M2  cotidal  chart.  In 
addition,  both  the  amplitudes  of  M2  and  Ki  are  large  near  the  outflow  of  the  Yangtze  and  Yellow 
Rivers  on  the  Chinese  coast.  Qualitative  comparisons  of  the  cotidal  charts  in  Figs.  2  and  3,  derived 
from  CU  model  results  for  M2  and  Ki,  with  previously  published  cotidal  charts  in  the  Yellow  Sea 
(e.g.,  Kang  et  al.  1990;  Choi  1984),  indicate  good  agreement. 


Table  1  —  Mean  Amplitudes  (m)  for  Eight  Primary  Tidal  Constituents 


RUN 

M2 

S2 

N2 

K2 

K, 

Oi 

Pi 

Qi 

NA 

0.7804 

0.3589 

0.1981 

0.0912 

0.2188 

0.1943 

0.0775 

0.0453 

A 

0.7769 

0.3011 

0.1161 

0.1010 

0.2165 

0.1696 

0.0765 

0.0343 
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Fig.  2  —  Cotidal  chart  of  the  M2  frequency  in  the  Yellow  and  East  China  Seas. 
Amplitudes  shown  as  solid  lines  at  20-in  intervals  and  dotted  lines  indicate 
phases  every  30  deg. 


5.0  INFLUENCE  OF  DATA-ASSIMILATION  STATION  LOCATION 
5.1  Experiments 

Both  TOPEX/Poseidon  altimeter  and  tide  gauge  data  are  available  in  the  Yellow  and  East 
China  Seas  and  can  be  assimilated  into  the  CU  tidal  model.  A  series  of  experiments  has  been 
designed  to  test  the  influence  of  the  location  of  assimilated  data  on  predictions  of  the  dominant  tidal 
constituents  M2  and  Ki.  The  role  of  data  assimilation  in  the  CU  model  is  assessed  to  facilitate  setup 
of  the  model  in  other  regions  and  to  assist  in  the  evaluation  of  tidal  computations.  Ideally,  one 
would  like  to  obtain  reasonable  tidal  predictions  using  the  least  amount  of  information  a  priori. 


118°  120°  122°  124°  126°  128°  130°E 

Fig.  3  —  Cotidal  chart  of  the  Kj  frequency  in  the  Yellow  and  East  China  Seas. 
Amplitudes  shown  as  solid  lines  at  5-m  intervals  and  dotted  lines  indicate  phases 
every  30  deg. 


Tidal  data  located  at  114  stations  are  divided  into  four  groups  delineated  by  their  location’s 
bathymetry  and/or  proximity  to  the  open  boundaries,  shown  in  Fig.  4,  and  named  accordingly:  ^ 

coastal,  shallow,  deep,  and  boundary  stations.  Coastal  stations  are  those  stations  that  lie  along  the 
shoreline.  These  stations  are  most  likely  to  record  highly  localized,  nonlinear  tidal  phenomena. 

Shallow  stations  are  those  stations  located  offshore  in  water  of  depths  less  than  200  m.  Data 

stations  whose  water  depth  is  greater  than  200  m  are  categorized  as  deep  stations.  Observations  that 

lie  closest  to  the  southern,  eastern,  and  Tsushima  Strait  open-ocean  boundaries  fall  into  the  boundary 

station  group.  Note  that  each  of  the  stations  in  the  boundary  group  is  taken  from  one  of  the  three  # 

coastal,  shallow,  or  deep  station  groupings.  Table  2  summarizes  the  number  of  stations  and 

distinguishing  location  of  the  station  groupings. 
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Table  2  —  Characteristics  of  the  Location  of  Assimilated  Data 


STATION  GROUP  TYPE 

TOTAL  NUMBER  OF  STATIONS 

RELATIVE  LOCATION 

All 

114 

Coastal 

51 

shoreline  (2  m  or  less) 

Shallow 

38 

2  m  <  depth  <  200  m 

Deep 

25 

depth  >  200  m 

Boundary 

17 

near  open  boundaries 

40° 

N 


38° 


36° 


34° 


32° 


30° 


28° 


26° 


Fig.  4  —  Four  data  assimilation  station  groupings — coastal,  shallow,  deep, 
and  boundary — located  throughout  the  Yellow  and  East  China  Seas 
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Six  model  simulations  are  conducted  for  both  the  M2  and  Ki  tidal  constituents  with  assimilation 
of  either  none,  all,  coastal,  shallow,  deep,  or  boundary  stations.  Cotidal  charts  for  each  of  the  model 
experiments  are  constructed  and  differences  in  the  co-amplitudes  between  tidal  heights  obtained 
by  assimilating  each  of  the  station  groupings  are  examined  to  determine  if  one  location  is  more 
influential  in  the  computation  of  primary  tides  than  another. 


5.2  Discussion  of  Results 

The  contour  maps  of  Fig.  5b-d  show  deviations  between  the  amplitude  of  the  M2  tide  obtained 
by  assimilation  of  all  observations  (Fig.  5a)  and  amplitudes  computed  by  assimilating  deep,  shallow, 
or  coastal  data  stations.  Positive  values  (solid  lines)  indicate  an  overprediction  and  negative  values 
(dotted  lines)  an  underprediction,  with  respect  to  the  “all-station”  case  (Fig.  5a).  In  Fig.  5b,  the 
largest  discrepancies  in  predicted  tides  occur  in  the  mid-Yellow  Sea,  along  the  east  and  west  lateral 
boundaries,  at  the  northernmost  bay  in  the  Bohai  Sea,  and  in  the  Taiwan  Strait.  In  all  of  these  areas 
except  the  upper  Bohai  Sea  and  the  Taiwan  Strait,  underprediction  is  dominant  and  ranges  between 
20  to  60  cm.  Overprediction  at  the  southern  end  of  the  Yellow  Sea  results  from  a  southeast  shift 
in  the  amphidrome  position.  When  only  those  stations  designated  as  boundary  station  data  in  Fig.  4 
are  assimilated,  the  pattern  of  differences  closely  follows  that  of  the  deep  station  assimilation, 
except  overprediction  at  the  Taiwan  Strait  and  underprediction  at  the  Tsushima  Strait  are  somewhat 
reduced.  Noticeable  improvements  result  from  the  assimilation  of  a  single  observation  in  each  of 
these  straits. 

Tidal  heights  for  M2  computed  using  only  shallow  observations,  as  shown  in  Fig.  5c,  differ 
only  slightly  from  the  solution  obtained  when  all  data  stations  are  assimilated.  Small  deviations 
occur  near  the  Yangtze  River,  in  isolated  near-shore  regions  along  the  west  Korean  coast,  and  in 
the  very  shallow  upper  Bohai  Sea.  Shallow  observations  near  the  Taiwan  and  Tsushima  Straits  are 
included  in  the  data  assimilation  and,  thus,  account  for  the  minimal  differences  seen  at  these  open 
boundaries.  Small  differences  between  predicted  M2  tidal  elevations  using  all  stations  versus  only 
shallow  stations  indicate  that  assimilation  of  shallow  observations  dominate  computation  of  the  M2 
tidal  constituent. 

The  assimilation  of  only  coastal  station  data  into  computations  of  the  M2  tide  produces  large 
variations  of  the  M2  amplitude  in  the  upper  Yellow  Sea  seen  in  Fig.  5d.  In  particular,  excessive 
underprediction,  up  to  60  cm,  is  located  offshore  from  the  Chinese  border.  Within  the  Bohai  Sea 
itself  differences  are  relatively  small.  Coastal  stations  bounding  the  Taiwan  Strait  fail  to  constrain 
M2  elevations  in  that  region.  Generally,  the  M2  tide  computed  in  shallow  regions  away  from  the 
coastline  does  not  agree  with  elevations  obtained  when  all  stations  are  assimilated.  Coastal  stations 
may  not  be  as  effective  in  predicting  the  M2  tide  in  part  because  stations  are  located  in  the  very 
nearshore,  outside  the  domain  modeled  due  to  the  minimum  bathymetry  specification  of  5  m. 
Furthermore,  coastal  stations  capture  highly  nonlinear,  localized  behavior  near  inlets  and  harbors 
and  may  not  be  well-suited  to  the  prediction  of  primary  tides. 

Similarly,  for  the  Ki  tidal  constituent,  contour  maps  in  Fig.  6b-d  show  differences  in  tidal 
elevation  computed  using  each  of  the  data-assimilation  groups,  deep,  shallow,  and  coastal,  as 
compared  to  the  Ki  solution  obtained  from  an  assimilation  of  all  stations  (Fig.  6a).  Again,  overprediction 
is  indicated  as  a  solid  line  and  underprediction  by  a  dotted  line.  As  shown  for  the  M2  tide,  inde¬ 
pendent  assimilation  of  the  deep  and  boundary  station  data  yields  interchangeable  results.  The  most 
significant  differences  in  tidal  height  predictions  in  comparison  to  the  all-data-assimilation  case  are 
concentrated  in  the  far  upper  Yellow  Sea,  seen  in  Fig.  6b.  Magnitudes  of  the  differences  lie  in  the 
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118°  120°  122°  124°  126°  128°  130° E  118°  120°  122°  124°  126°  128°  130°  E 

Fig.  6  —  (a)  Cotidal  chart  of  the  Kj  tide  with  amplitudes  in  5-cm  solid  contours  and  phases  in  dotted  lines 
at  30-deg  intervals.  Differences  between  the  Kj  amplitude  computed  by  assimilating  all  data  stations  and 
only  (b)  deep  stations,  (c)  shallow  stations,  or  (d)  coastal  stations.  Contours  are  plotted  at  1-cm  intervals 
with  overprediction  (positive)  as  solid  lines  and  underprediction  (negative)  as  dotted  lines. 
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range  of  8  cm  or  less  with  exceptions  in  the  northern  Bohai  sea,  at  the  Yangtze  River  entrance,  at 
Incheon  on  the  South  Korean  coast,  and  in  the  Taiwan  Strait  where  differences  reach  up  to  20  cm. 

The  assimilation  of  shallow  observations  results  in  a  solution  for  Ki  that  closely  mirrors  that 
obtained  by  assimilation  of  all  observations  as  shown  in  Fig.  6c.  Locations  near  Incheon  and  the 
Yangtze  River  offer  the  only  major  discrepancies.  Again,  the  assimilation  of  shallow  data  dominates 
the  computed  tidal  solution,  in  this  case  for  the  Ki  tide. 

The  assimilation  of  coastal  station  data  leads  to  a  distinctive  over-  and  underprediction  north- 
south  division  in  the  Yellow  Sea  as  depicted  in  the  difference  map  of  Fig.  6d.  Underprediction  is 
indicated  along  the  west  Korean  coast  and  overprediction  along  the  Chinese  coast.  This  pronounced 
feature  can  be  attributed  to  a  shift  in  the  centrally  located  amphidrome  in  the  Ki  tide.  Magnitudes 
of  the  differences  remain  less  than  10  cm  except  at  points  near  the  Taiwan  Strait,  on  the  Chinese 
coast,  an  in  the  extremely  shallow  regions  of  the  Bohai  Sea. 

Tidal  predictions  for  Ki,  which  independently  assimilate  the  four  designated  groups  of  data 
stations,  exhibit  only  small  changes  between  them  as  compared  to  the  larger  differences  associated 
with  computations  of  the  M2  tide  using  each  of  the  station  groupings. 


5.3  Conclusions 

The  assimilation  of  data  into  the  CU  tidal  model  is  influential  in  predicting  the  tidal  response 
of  the  Yellow  Sea  basin.  From  an  analysis  of  four  data  station  groupings,  coastal,  shallow,  deep, 
and  open  boundary  located  stations,  observations  located  on  the  continental  shelf  away  from  the 
coastline  are  most  influential  in  the  prediction  of  primary  tides,  M2  and  Kj,  in  the  Yellow  Sea. 
Effects  of  deep-water  stations  or  those  located  near  large,  open-ocean  boundaries  do  not  reach 
beyond  the  deeper  waters  and  have  little  impact  on  tidal  computations  on  the  continental  shelf. 
However,  the  analysis  here  has  shown  that  tidal  data  assimilated  in  the  open  water  of  narrow  straits 
can  improve  local  tidal  height  computations.  The  influence  of  coastal  station  data  was  largest  in 
small  near-shore  regions  and  was  not  a  major  factor  in  the  computation  of  the  primary  tides,  M2 
and  Ki-  Coastal  gauges  are  often  located  in  harbors  or  inlets  where  nonlinear  tidal  interactions  are 
prominent  and  wave  action  may  be  present,  a  phenomena  not  included  in  the  model  physics.  A 
preliminary  examination  of  the  M4  overtide  solution  (not  shown)  indicates  that  assimilation  of 
coastal  observations  may  be  more  significant  in  predicting  nonlinear  tidal  patterns  in  the  northern 
Yellow  Sea.  A  thorough  look  at  the  role  of  data  assimilation  in  the  prediction  of  nonlinear  tidal 
constituents  is  slated  for  future  study. 

This  documentation  of  CU  model  performance  should  enhance  development  of  data  sampling 
strategies  to  improve  future  tidal  model  predictions.  Issues  of  how  often  to  assimilate  data  or  how 
to  assign  appropriate  weights  to  each  observation  are  not  addressed  here. 


6.0  TIDAL  MODEL  EVALUATION 

From  the  original  114  data  locations  in  the  Yellow  Sea  (Fig.  4),  49  points  are  selected  to  serve 
as  measurements  against  which  tidal  model  computations  are  evaluated.  These  observation  locations, 
depicted  in  Fig.  7,  are  selected  so  that  their  spatial  distribution  covers  much  of  the  Yellow  and  East 
China  Seas  model  domain.  At  the  remaining  65  locations,  also  shown  in  Fig.  7,  measured  tidal 
heights  are  assimilated  into  tidal  model  computations  at  each  time  step,  intervals  of  24  s. 


14 


Cheryl  Ann  Blain 


Tidal  elevations  are  assimilated  as  amplitudes  and  phases  of  each  tidal  constituent.  Weights  assigned 

to  the  stations  are  fixed  at  90%.  ^ 

The  CU  tidal  model  is  evaluated  by  comparing  computations  of  the  amplitudes  and  phases  of 
eight  major  tidal  constituents,  M2,  S2,  N2,  K2  and  Ki,  Oi,  Pi,  Qi,  to  measured  tidal  amplitudes  and 
phases  at  the  49  measurement  locations  throughout  the  Yellow  and  East  China  Seas.  For  compo¬ 
nents  N2,  K2,  Pi,  and  Qi,  observations  are  only  available  at  20,  17,  21,  and  12  of  the  49  measurement 
locations,  respectively.  Two  simulations.  Run  A  and  Run  NA  mentioned  previously,  are  conducted  # 

for  this  model  evaluation.  Recall  Run  A  includes  the  assimilation  of  observed  tidal  heights  at 
65  stations,  whereas  Run  NA  has  no  data  assimilation. 
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6.1  Comparison  to  Observations 

From  Foreman  (1993),  differences  between  model-derived  and  measured  tidal  amplitudes  and 
phases  are  computed  as  distances,  D,  in  the  complex  plane; 

D  =  {(Ao  cosPo  -  Am  cos  Pm)^  +  (A^  sinPo  -  A^  sin  Pm)^}^'^,  (8) 

where  Ao,  Am,  Po,  and  Pm  are  the  observed  and  modeled  amplitudes  and  phases,  respectively. 
Average  values  of  D  for  each  tidal  constituent,  as  well  as  average  differences  between  model 
computed  and  observed  amplitudes  and  phases,  are  given  in  Table  3  for  Run  A,  the  data-assimilative 
simulation.  From  Table  3,  the  average  M2  amplitude  error  is  approximately  13  cm  and  the  distance 
in  the  complex  plane  is  nearly  30  cm,  a  significant  portion  of  which  is  due  to  large  phase  errors. 
Though  the  largest  discrepancies  between  model  results  and  observations  occur  in  the  prediction  of 
the  M2  tide,  as  measured  by  D,  proportionally,  the  differences  in  amplitude  for  M2  and  Ki  are  less 
than  13%  of  mean  measured  amplitudes  for  these  constituents.  Model  performance  is  considerably 
worse  for  the  four  less  prominent  tidal  frequencies,  N2,  K2,  Pi,  and  Qi,  and  may  be  due  in  part  to 
the  reduced  number  of  measurement  points  available  for  comparison. 

Spatially,  the  most  significant  differences  between  model  computed  and  observed  tidal  amplitude 
and  phase  are  concentrated  along  the  west  coast  of  Korea  and  in  the  Taiwan  Strait.  Near  Incheon, 
located  on  the  northwest  South  Korean  coast,  tide  model  results  are  extremely  poor  in  comparison 
to  observed  values.  The  complex  coastline  and  very  shallow  depths  at  Incheon  are  not  well  repre¬ 
sented  by  the  coarse  CU  Yellow  Sea  model  and  largely  account  for  the  significant  disparity.  Some 
errors  in  the  shallowest  water  may  be  attributed  to  coastal  gauges  located  in  harbors  or  inlets  where 
nonlinear  tidal  interactions  are  prominent  and  wave  action  may  be  present,  a  phenomena  not  included 
in  the  model  physics.  Furthermore,  coastal  data  is  assimilated  into  the  model  at  nearest  nodal 
locations  introducing  an  additional  source  of  error. 

In  the  deep  waters,  particularly  those  locations  where  TOPEX/Poseidon  altimetry  data  has  been 
assimilated,  model  predictions  are  quite  good.  Tidal  data  assimilated  in  open  waters  of  narrow 
straits  improves  local  tidal  height  computations.  Deep-water  assimilation  of  data  has  little  impact 
on  tidal  predictions  on  the  continental  shelf. 


Table  3  —  Absolute  Errors  from  Run  A  for  Eight  Primary  Tidal  Constituents 


ABSOLUTE 

ERRORS 

M2 

S2 

N2 

K2 

Ki 

Oi 

Pi 

Qi 

D,  m 

0.290 

0.153 

0.101 

0.083 

0.064 

0.052 

0.033 

0.027 

\Ao  -  Am\,  m 

0.126 

0.071 

0.061 

0.066 

0.031 

0.028 

0.026 

0.16 

\P 0  ~  P ml.  deg 

46.18 

82.16 

38.66 

76.02 

26.48 

12.11 

28.01 

31.37 

Measured  Mean 
Amplitude 

0.9880 

0.3471 

0.1836 

0.1134 

0.2606 

0.1955 

0.0773 

0.0336 
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Finally,  RMS  errors  are  computed  for  each  major  constituent  using  the  standard  formula: 


where  i  is  the  tidal  constituent,  tm,  to  are  the  model  and  observed  amplitudes  or  phases,  and  n  is 
the  number  of  observation  points.  RMS  errors  computed  for  Runs  NA  and  A  are  presented  in 
Table  4a  for  amplitude  and  in  Table  4b  for  phase.  A  one-  to  two-fold  reduction  in  the  RMS  error 
is  evident  in  tidal  amplitudes  computed  using  a  data-assimilative  model,  i.e..  Run  A,  as  compared 
to  the  no  assimilation  case.  Run  NA. 

RMS  amplitude  errors  represented  as  a  percentage  of  the  mean  measured  amplitude  of  each 
constituent  are  also  given  in  Table  4a.  For  Run  A,  the  lowest  percentage  RMS  errors  for  amplitude, 
17-19%,  are  linked  to  the  diurnal  tidal  constituents  Ki  and  Oi.  Semidiurnal  components  M2  and 
S2  have  proportional  RMS  amplitude  errors  ranging  from  30-40%.  Significantly  larger  RMS 
percentage  amplitude  errors  are  computed  for  remaining  constituents  and  may  again  be  a  consequence 
of  fewer  observational  points  available  for  comparison.  RMS  phase  errors  are  noticeably  poor  with 
an  RMS  phase  difference  just  over  110°  for  the  M2  tide  and  nearly  87°  for  the  Ki  tide. 


6.2  Sensitivity  of  Bottom  Friction  Coefficient 

Table  5  shows  the  effect  of  varying  c/  on  mean  tidal  amplitudes,  as  well  as  RMS  and  absolute 
errors  of  amplitude  and  phase,  for  the  major  semidiurnal  and  diurnal  frequencies,  M2  and  Ki. 
Selected  friction  factors  in  the  range  from  0.0005  to  0.004  are  fixed  over  the  entire  Yellow  and  East 
China  Seas  model  domain.  An  increase  in  the  friction  factor  most  notably  affects  model  computations 
of  amplitude  and  phase  for  the  M2  tidal  constituent.  The  mean  M2  amplitudes  and  RMS  amplitude 
errors  decrease  with  increasing  c/ values  and  are  reduced  by  approximately  10  and  12%,  respectively, 
over  the  range  of  friction  factors  considered.  RMS  phase  errors  computed  for  each  friction  coefficient 
are  within  10°  with  a  minimum  occurring  at  c/=  0.002. 


Table  4a  —  RMS  Amplitude  Errors  (m)  for  Eight  Primary  Tidal  Constituents 


RUN 

M2 

S2 

N2 

K2 

Ki 

Oi 

Pi 

Qi 

NA 

0.4493 

0.2604 

0.1440 

0.1133 

0.0748 

0.0788 

0.0519 

0.0516 

A 

0.2989 

0.1331 

0.0996 

0.1093 

0.0439 

0.0363 

0.0551 

0.0285 

%  of  Obs. 

30 

38 

54 

96 

17 

19 

71 

85 

Table  4b  —  RMS  Phase  Errors  (deg)  for  Eight  Primary  Tidal 

Constituents 


RUN 

M2 

S2 

N2 

K2 

Ki 

Oi 

Pt 

Qi 

NA 

136.33 

141.21 

83.00 

137.02 

71.90 

71.37 

99.06 

54.20 

A 

110.31 

151.82 

83.54 

142.50 

86.91 

49.35 

105.43 

41.89 
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Table  5a  —  Run  A  Statistics  for  the  M2  Tide  Under  a  Varying  Friction  Factor 


FRICTION 

MEAN 

RMS  ERRORS 

IAq  Affi\ 

\Po  -  Pm\ 

FACTOR,  Cf 

AMP,  m 

AMP,  m 

PHASE,  deg 

Z),  m 

m 

deg 

0.0005 

0.8282 

0.2916 

121.04 

0.291 

0.169 

52.33 

0.001 

0.8056 

0.2515 

121.93 

0.279 

0.148 

52.74 

0.002 

0.7769 

0.2040 

110.59 

0.290 

0.126 

46.18 

0.003 

0.7574 

0.1851 

111.09 

0.285 

0.115 

46.19 

0.004 

0.7418 

0.1817 

112.28 

0.287 

0.116 

46.89 

Table  5b  —  Run  A  Statistics  for  the  Kj  Tide  Under  a  Varying  Friction  Factor 


FRICTION 

MEAN 

RMS  ERRORS 

\Aq  Affil 

\Po  -  Pm\ 

FACTOR,  Cf 

AMP,  m 

AMP,  m 

PHASE,  deg 

D,  m 

m 

deg 

0.0005 

0.2162 

0.0622 

72.33 

0.080 

0.040 

26.86 

0.001 

0.2159 

0.0463 

75.12 

0.069 

0.034 

26.94 

0.002 

0.2165 

0.0435 

76.17 

0.064 

0.031 

26.48 

0.003 

0.2165 

0.0428 

76.82 

0.062 

0.031 

26.53 

0.004 

0.2162 

0.0423 

77.05 

0.060 

0.031 

26.46 

The  varying  bottom  friction  coefficient  effects  the  Ki  frequency  to  a  much  lesser  extent.  Ki 
mean  amplitudes  remain  close  in  magnitude,  changing  less  than  0.2%,  and  exhibit  no  clear  trend 
in  relation  to  the  friction  factor  specification.  In  contrast,  RMS  amplitude  errors  computed  for  the 
Ki  tide  decrease  with  increasing  Cf  values  and  RMS  phase  errors  exhibit  the  reverse.  Keep  in  mind 
that  these  deviations  in  model  computations  of  the  K]  tide  in  response  to  changing  friction 
coefficients  are  slight.  The  ratio  of  M2  to  Kj  amplitudes  is  relatively  stable,  changing  by  only  3%, 
regardless  of  the  bottom  friction  factor  selected. 

Absolute  error  measures  in  Table  5  echo  the  M2  RMS  amplitude  error  trend  that  differences 
decrease  for  increasing  values  of  Cf.  Again,  changes  in  errors  associated  with  the  M2  tide  are  more 
dramatic  than  those  of  the  Ki  tide.  Amplitude  and  phase  absolute  error  measures  for  M2  and  Ki 
tidal  frequencies,  for  the  most  part,  reach  a  minimum  when  the  bottom  friction  factor  is  in  the  range 
of  0.002-0.003.  These  c/  values  correspond  to  cf  =  0.0025,  a  value  most  often  reported  in  the 
literature  for  Yellow  Sea  tidal  simulations  (i.e.,  Kang  et  al.  1990;  Lee  and  Beardsley  1993;  Zhenxing 
et  al.  1994). 


6.3  Conclusions 

An  evaluation  of  the  CU  tidal  model  computations  for  eight  primary  tidal  constituents  demonstrates 
that  the  assimilation  of  data  into  tidal  models  can  be  influential  in  predicting  the  tidal  response  of 
a  basin.  A  coarse  tidal  model  of  the  Yellow  Sea  produces  reasonable  computations  of  tidal  elevation 
with  low  RMS  errors  when  compared  at  49  observation  sites.  RMS  errors  derived  from  CU  model 


simulations  of  tidal  elevation  computed  with  and  without  data  assimilation  indicate  that  assimilation 
of  tidal  information  reduces  error  by  a  factor  of  two.  An  abundance  of  reliable  sea-surface  height 
observations  from  coastal  stations  and  satellite  sources  in  the  Yellow  Sea  partially  compensates  for 
the  coarse  grid  resolution  of  the  CU  tidal  model. 

The  M2  tidal  constituent  is  the  dominant  frequency  in  the  Yellow  Sea.  Amplitudes  reach  nearly 
2  m  along  the  coast  of  west  Korea.  Ki  is  the  most  significant  diurnal  component  with  amplitudes 
just  over  1/4  of  those  for  the  M2  tide.  Computed  amplitudes  of  the  M2  and  Ki  tides  are  within  13% 
of  measured  tidal  elevations  at  49  locations.  Phases  computed  by  the  CU  Yellow  Sea  tidal  model 
are  poorly  represented  and  warrant  further  investigation.  Additional  improvements  would  no  doubt 
result  from  using  a  more  highly  resolved  grid  discretization. 

A  sensitivity  study  of  bottom  friction  coefficient  specification  on  model  results  suggests  the 
friction  factor  is  most  influential  in  computations  of  the  M2  tide.  As  expected,  an  increasing  friction 
coefficient  decreases  tidal  amplitudes.  Absolute  errors  reveal  an  optimal  friction  coefficient  ranging 
between  0.002  and  0.003,  values  that  agree  well  with  published  numbers  for  the  Yellow  Sea. 

Generally,  tidal  amplitudes  for  all  eight  constituents  computed  using  data  assimilation  are 
overdamped.  This  overdamping  of  the  tidal  amplitude  may  be  the  result  of  numerical  deficiencies 
in  the  model  or  inappropriate  specification  of  the  horizontal  mixing  parameter,  since  the  bottom 
friction  coefficient  has  been  shown  to  have  limited  influence.  Issues  of  how  often  to  assimilate  data 
or  how  to  assign  appropriate  weights  to  each  observation  are  not  addressed  here.  Nonlinear  tidal 
interactions  have  been  ignored  in  this  paper  but  are  the  focus  of  current  studies  using  the  CU 
Yellow  Sea  tidal  model. 


7.0  SUMMARY 

The  Yellow  Sea  is  a  very  shallow-water  body  with  a  strong  tidal  signal.  Its  strategic  importance 
had  led  the  U.  S.  Navy  to  study  tidal  response  in  the  Yellow  and  East  China  Seas  using  the 
Colorado  University  data-assimilative,  barotropic  tidal  model.  The  CU  model  has  been  adapted  to 
accept  open-ocean  boundary  forcing  from  the  Grenoble  global  tidal  data  base,  which  is  considered 
more  reliable  than  tides  currently  obtained  from  the  CU  global  tidal  model.  Furthermore,  the 
number  and  quality  of  the  data-assimilation  stations  utilized  by  the  CU  Yellow  Sea  model  is 
improved.  A  detailed,  semi-automated  procedure  for  assigning  data  stations  for  assimilation  into  the 
CU  model  domain  has  been  developed. 

Characterization  of  the  tidal  elevations  in  the  Yellow  and  East  China  Seas  from  CU  model 
simulations  demonstrate  that  the  M2  constituent  is  the  dominant  tidal  frequency  in  the  Yellow  Sea 
with  amplitudes  of  nearly  2  m  along  the  coast  of  west  Korea.  Ki  is  the  most  significant  diurnal 
component  with  amplitudes  1/4  of  the  M2  tide. 

A  series  of  experiments  is  conducted  to  test  the  influence  of  the  location  of  assimilated  data 
on  predictions  of  the  dominate  tidal  constituents,  M2  and  Kj.  From  an  analysis  of  four  data  station 
groupings,  coastal,  shallow,  deep,  and  open  boundary  located  stations,  observations  located  on  the 
continental  shelf  away  from  the  coastline,  are  most  influential  in  the  prediction  of  primary  tides, 
M2  and  Ki,  in  the  Yellow  Sea.  Coastal  station  data  may  be  more  influential  in  the  prediction  of 
nonlinear  tidal  interactions  since  these  gauges  are  often  located  in  harbors  or  inlets  where  nonlinear 
tidal  interactions  are  prominent  and  wave  action  may  be  present.  Furthermore,  coastal  data  is 
assimilated  into  the  model  at  the  nearest  nodal  locations  introducing  an  additional  source  of  error. 
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This  examination  of  the  data-assimilation  component  of  the  CU  model  should  enhance  development 
of  data  sampling  strategies  to  improve  future  tidal  model  predictions. 

Validation  of  the  CU  tidal  model  is  accomplished  by  comparing  predicted  and  observed  tidal 
elevations  at  49  locations.  For  eight  primary  tidal  constituents,  absolute  and  RMS  amplitude  and 
phase  errors  are  examined  for  data-assimilative  and  nonassimilative  computations.  CU  model  results 
for  the  M2  and  Ki  tides  are  within  13%  of  measured  values,  but  phases  are  poorly  represented. 
RMS  errors  for  the  data-assimilative  simulation  are  lower  than  those  for  the  nonassimilative  run, 
but  comparison  of  absolute  tidal  amplitude  errors  indicate  overdamping  by  the  CU  model. 

A  look  at  the  sensitivity  of  the  bottom  friction  coefficient  specification  on  model  computations 
shows  that  the  M2  tide  is  most  affected.  An  increasing  friction  coefficient  decreases  tidal  amplitudes 
although  absolute  errors  indicate  an  optimal  friction  factor  between  0.002  and  0.003,  which  is  in 
good  agreement  with  published  values  for  the  Yellow  Sea. 
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